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ABSTRACT 


A procedure  has  been  developed,  based  on  a dilute-solution  model  of  the  liquid  phase  and  the 
rigorous  virial  expansion  for  the  vapor  phase,  for  modeling  the  vapor-liquid  equilibrium  between 
water  and  a synthesis  gas  consisting  primarily  of  H2,  CO,  and  CO2.  Calculations  of  the  vapor- 
liquid  equilibrium  have  been  performed  at  temperatures  and  pressures  typical  for  a quench  step 
in  an  IGCC  power  plant.  The  uncertainty  in  the  calculated  vapor-phase  mole  fraction  of  water  is 
less  than  0.01,  and  is  due  primarily  to  the  omission  of  higher-order  terms  in  the  virial  expansion. 

Other  thermodynamic  calculation  methods  were  examined  and  compared  to  the  more  rigorous 
calculations.  The  ideal-gas  assumption  for  the  vapor  phase  seriously  underestimates  (by  about 
0.09  mole  fraction  at  typical  conditions)  the  water  content  of  the  equilibrium  vapor.  The  Peng- 
Robinson  equation  performs  much  better,  producing  vapor-phase  mole  fractions  of  water  that 
differ  from  the  more  rigorous  results  by  amounts  similar  to  the  uncertainty  in  the  virial  approach. 

This  work  was  performed  for  the  Program  for  Technology  Innovation  of  the  Electric  Power 
Research  Institute  (EPRI).  A version  of  this  report  with  additional  process-related  comments 
from  the  sponsor  perspective  is  being  published  separately  as  an  EPRI  Report. 
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1 

INTRODUCTION 


Background 

The  purpose  of  this  work  is  to  provide  data  relevant  to  the  design  of  Integrated  Gasification 
Combined  Cycle  (IGCC)  power  plants.  In  some  proposed  designs,  a hot  synthesis  gas,  consisting 
mostly  of  H2,  CO,  and  CO2,  is  quenched  with  water.  The  amount  of  water  in  the  resulting  gas  is  a 
significant  factor  in  the  economics  of  the  cycle.  Previous  efforts  to  model  this  step,  using  a 
variety  of  thermodynamic  models,  have  come  up  with  significantly  different  results.  This  project 
was  undertaken  to  perform  calculations  with  a rigorous  thermodynamic  model  in  order  to 
provide  reliable  data  for  the  water  vapor  content,  with  reasonable  estimates  of  uncertainties. 


Problem  Definition 

The  calculations  performed  consist  of  the  computation  of  thermodynamic  equilibrium  between 
vapor  and  liquid  (vapor-liquid  equilibrium,  or  VLE).  This  equilibration  takes  place  at  a specified 
temperature  T and  absolute  pressure  p.  Based  on  input  from  the  sponsor,  a “baseline”  case  was 
defined  at  a temperature  of  495  K (approximately  431  °F)  and  a pressure  of  4.5  MPa 
(approximately  650  psia). 

The  “baseline”  composition  for  the  dry  gas  was  also  taken  from  sponsor  input  and  is  shown  in 
Table  1-1.  Other,  trace  components  in  the  gas  phase  may  exist  (such  as  H2S,  NH3,  and  CH4),  but 
they  are  outside  the  scope  of  this  project.  While  the  main  result  of  interest  (mole  fraction  of 
water  in  equilibrium  vapor)  is  not  very  sensitive  to  the  total  amount  of  water  in  the  system,  for 
consistency  all  calculations  assumed  five  times  as  many  moles  of  water  as  moles  of  dry  gas. 


Table  1-1 

Baseline  Composition  of  Synthesis  Gas 


Gas 

Mole  Fraction 

Ar 

0.01 

H2 

0.39 

N2 

0.01 

CO 

0.41 

0 

0 

ro 

0.18 
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Approach 

In  order  to  solve  the  VLE  problem,  the  fugacity  of  each  component  in  the  liquid  phase  is  equated 
with  its  fugacity  in  the  vapor  phase  under  the  constraint  of  fixed  T and  p.  In  the  terminology  of 
chemical  engineering,  this  is  known  as  a “flash”  calculation.  The  computational  algorithm  for 
solving  the  flash  problem  is  described  in  Appendix  A.  This  algorithm  requires  a model  for  the 
liquid-phase  and  vapor-phase  fugacities  of  each  component  as  a function  of  temperature, 
pressure,  and  composition. 

For  modeling  the  fugacities  in  the  liquid  phase  at  these  moderate  pressures,  where  the  solubility 
of  the  gases  in  water  is  small,  it  is  sufficient  to  invoke  the  dilute-solution  approximation  known 
as  Henry’s  Law.  This  primarily  requires  knowing  the  Henry’s  constant  for  each  gas  in  water  as  a 
function  of  temperature.  Additional  data  on  the  partial  molar  volumes  of  gases  in  solution 
provide  a small  correction  for  the  effect  of  pressure  on  the  fugacity.  The  liquid-phase  fugacity 
model  used  in  this  work  is  completely  documented  in  Appendix  B. 

The  most  important  factor  in  determining  the  water  content  of  the  equilibrium  vapor  is  the  vapor- 
phase  fugacity  model.  At  the  conditions  of  interest,  an  ideal-gas  assumption  (common  in  many 
engineering  calculations)  would  seriously  overestimate  the  fugacity  of  water,  leading  to  an 
underestimation  of  the  equilibrium  water  mole  fraction.  The  equation-of-state  approaches 
commonly  used  in  engineering  thermodynamics  are  designed  primarily  for  hydrocarbon 
mixtures,  so  they  cannot  be  expected  to  be  reliable  for  mixtures  of  water  with  nonpolar  gases  as 
encountered  in  this  work. 

Our  approach  here  is  to  model  the  vapor-phase  fugacities  by  the  virial  expansion,  which  is  a 
rigorous  series  of  corrections  to  the  ideal-gas  law.  At  the  moderate  pressures  of  interest  here,  the 
second  virial  coefficient  (first  correction  to  the  ideal  gas)  provides  sufficient  accuracy.  Therefore, 
if  reliable  values  of  the  second  virial  coefficients  among  all  components  in  the  mixture  can  be 
obtained,  accurate  values  of  the  vapor-phase  fugacities  can  be  obtained.  Appendix  C documents 
the  second-virial-coefficient  data  used  in  this  work,  and  their  use  in  fugacity  calculations. 

For  purposes  of  comparison,  we  also  performed  vapor-phase  fugacity  calculations  with  the  ideal- 
gas  model,  and  with  a common  engineering  equation  of  state  model  (the  Peng-Robinson 
equation,  applied  only  to  the  vapor  phase).  These  models  are  also  described  in  Appendix  C. 
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RESULTS 


Baseline  Case  Results 

Table  2-1  shows  the  results  of  VLE  calculations  for  the  baseline  case  as  defined  in  Chapter  1 
(495  K,  4.5  MPa),  for  all  three  of  the  vapor  fugacity  calculation  methods  considered  in  this  work. 
While  the  quantity  of  direct  interest  in  this  project  is  the  vapor-phase  mole  fraction  of  water,  all 
compositions  in  both  phases  are  shown  for  completeness. 

The  rigorous  calculation  from  virial  coefficients  produces  a vapor  mole  fraction  of  water  of 
0.5731,  much  larger  than  would  be  obtained  from  an  ideal-gas  calculation.  Perhaps  surprisingly, 
use  of  the  Peng-Robinson  equation  for  vapor  fugacities  produces  results  in  fairly  good  agreement 
with  the  more  rigorous  calculation. 

Table  2-1 

Equilibrium  Liquid  (x,)  and  Vapor  (y,)  Moie  Fractions  Calcuiated  for  the  Baseline  Case,  with 
Vapor  Fugacities  Computed  by  the  ideal-Gas  (IG),  Virial  (Vir),  and  Peng-Robinson  (PR) 
Methods 


/■ 

xm 

xXVir) 

x/(PR) 

ym 

y/(Vir) 

y/(PR) 

H20 

0.9986 

0.9988 

0.9989 

0.4817 

0.5731 

0.5666 

Ar 

6.87x1 

6.15x10“® 

5.96x10® 

0.0052 

0.0043 

0.0043 

H2 

2.98x10"* 

2.72x10^ 

2.68x10"’ 

0.2026 

0.1668 

0.1693 

N2 

4.89x10® 

4.38x10“® 

4.32x10  ® 

0.0052 

0.0043 

0.0043 

CO 

2.87x10“^ 

2.56x10"’ 

2.53x10^ 

0.2131 

0.1754 

0.1780 

CO2 

7.58x10“'’ 

6.32x10^ 

6.14x10^ 

0.0922 

0.0762 

0.0774 

Effects  of  Temperature  and  Pressure 

To  assess  the  effect  of  temperature  and  pressure  on  the  results,  calculations  were  performed  both 
over  a range  of  temperatures  (480  K to  510  K,  fixing  the  pressure  at  its  baseline  value  of  4.5 
MPa),  and  over  a range  of  pressures  (4.0  MPa  to  5.0  MPa,  fixing  the  temperature  at  its  baseline 
value  of  495  K).  The  results  of  these  calculations  are  shown  in  Figures  2-1  and  2-2. 
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Temperature  (K) 


Figure  2-1 

Effect  of  Temperature  on  Equilibrium  Mole  Fraction  of  Water  in  Vapor  at  4.5  MPa 


Pressure  (MPa) 
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Figure  2-2 

Effect  of  Pressure  on  Equilibrium  Mole  Fraction  of  Water  in  Vapor  at  495  K 
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Figure  2-1  shows  that,  as  would  be  expected  (because  of  the  increase  of  the  vapor  pressure  of 
water  with  temperature),  increasing  the  temperature  at  constant  pressure  results  in  more  water  in 
the  vapor  phase.  The  difference  between  the  virial  result  and  the  other  two  vapor  fugacity 
methods  grows  slowly  with  temperature,  but  over  the  entire  range  shown  the  ideal-gas  results  are 
seriously  in  error,  while  the  Peng-Robinson  results  are  reasonable. 

Figure  2-2  (note  that  the  scale  is  different  from  that  in  Figure  2-1)  shows  that  an  increase  in 
pressure  at  constant  temperature  results  in  a reduced  mole  fraction  of  water  in  the  vapor.  This 
also  makes  physical  sense;  the  partial  pressure  of  water  is  determined  primarily  by  its  vapor 
pressure  as  a function  of  temperature,  so  any  increase  in  pressure  must  come  from  increasing  the 
amount  of  other  components  in  the  synthesis  gas.  Over  this  range  of  pressures  at  fixed 
temperature,  the  relative  differences  among  the  three  vapor  fugacity  methods  is  nearly  constant. 


Estimation  of  Uncertainties 

Three  main  sources  of  uncertainty  in  the  calculated  vapor-phase  water  content  can  be  identified 
in  our  model.  In  the  discussion  below,  we  will  use  the  term  “uncertainty”  to  refer  to  an  expanded 
uncertainty  with  coverage  factor  k=2,  which  approximately  corresponds  to  a 95  % confidence 
interval. 

Uncertainty  from  Liquid  Fugacity  Modei 

The  liquid-phase  fugacity  model  introduces  only  negligible  uncertainty.  The  fugacity  of  pure 
water,  which  is  the  main  part  of  this  model  that  determines  the  water  content  in  the  vapor,  is 
determined  mainly  by  the  vapor  pressure  of  water,  which  has  an  uncertainty  of  approximately 
0.02  % in  this  temperature  range.  [1]  A small  correction  for  pressure  is  on  the  order  of  1 %,  and 
it  is  determined  by  the  liquid  molar  volume,  which  is  known  to  within  0.1  %.  [1] 

The  gas  solubility  in  the  liquid  phase  is  determined  primarily  by  the  Henry’s  constants  for  the 
different  gases  in  water,  which  have  uncertainties  of  roughly  5 % in  this  temperature  range.  [2] 
However,  because  the  resulting  solubilities  are  still  quite  small  (see  Table  2-1),  the  only 
significant  effect  of  this  uncertainty  on  the  VLE  results  is  on  the  mole  fractions  of  the  gases  in 
the  liquid  phase;  the  water  content  in  the  vapor  is  not  significantly  affected. 

Uncertainty  from  Second  Viriai  Coefficients 

The  water-gas  and  water-water  second  virial  coefficients  all  have  some  uncertainty  associated 
with  them  (see  references  in  Appendix  C).  We  tested  the  impact  of  this  uncertainty  on  the  results 
by  perturbing  each  of  these  virial  coefficients  (excluding  those  for  the  trace  gas  components 
argon  and  nitrogen)  by  its  (k=2)  uncertainty  and  repeating  the  calculations  for  the  baseline  case. 
The  results  are  shown  in  Table  2-2.  The  results  in  the  table  show  that  this  contribution  to  the 
uncertainty  is  also  relatively  small,  contributing  to  the  uncertainty  of  the  water  vapor  mole 
fraction  by  about  0.002  for  the  second  virial  coefficient  of  pure  water  and  by  less  than  0.001  for 
each  of  the  water-gas  coefficients. 
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Table  2-2 

Effect  on  Water  Content  in  Vapor  of  Perturbing  Second  Virial  Coefficient  Upward  by  its 
(lc=2)  Uncertainty 


Second  Virial  Coefficient  Perturbed 

Water  Mole  Fraction  in  Vapor 

(Baseline  case) 

0.5731 

H2O-H2O 

0.5708 

H2-H2O 

0.5729 

CO-H2O 

0.5724 

0 

CM 

1 

CM 

0 

0 

0.5728 

There  are  also  uncertainties  in  the  second  virial  coefficients  not  involving  water,  but  their 
influence  on  the  water  content  of  the  vapor  is  relatively  negligible. 

Uncertainty  from  Higher  Virial  Coefficients 

The  truncation  of  the  virial  expansion  after  the  first  correction  necessarily  introduces  some  error 
due  to  the  omission  of  higher-order  terms.  These  errors  will  increase  with  increasing  density.  We 
can  estimate  the  effect  of  higher-order  terms  by  considering  the  binary  system  H2O-CO2.  For  this 
system,  experimental  values  of  third  virial  coefficients  in  the  temperature  range  of  interest  are 
available  for  pure  water,  pure  CO2,  and  for  the  binary  interactions.  [1,3,4]  These  values  have 
large  uncertainties,  but  they  serve  the  purpose  of  estimating  the  magnitude  of  the  uncertainty  due 
to  virial  coefficients  beyond  the  second. 

We  establish  a “baseline”  case  for  the  H2O-CO2  mixture  at  the  same  T and  p as  our  main  case. 
For  this  case,  calculation  with  our  second-virial  model  yields  a water  vapor  mole  fraction  of 
0.5788.  Repeating  the  calculation  with  the  inclusion  of  third  virial  coefficients  produces  a vapor 
mole  fraction  of  0.5857.  Further  analysis  reveals  that  almost  all  of  this  change  is  due  to  the  third 
virial  coefficient  of  pure  water. 

We  therefore  conclude  that  the  omission  of  higher  virial  coefficients  causes  the  vapor  mole 
fraction  of  water  to  be  underestimated  by  an  amount  on  the  order  of  0.005  to  0.01.  This  is  the 
largest  contribution  to  the  uncertainty  in  our  calculation.  The  fact  that  most  of  this  deviation 
appears  to  be  due  to  the  third  virial  coefficient  of  pure  water  suggests  that  including  only  that 
coefficient  could  reduce  this  source  of  uncertainty.  However,  the  improvement  in  uncertainty 
would  be  limited  by  the  large  uncertainty  in  water’s  third  virial  coefficient;  two  sources  for  this 
coefficient  differ  from  each  other  by  about  25  %.  [1,5] 


Effect  of  Gas  Composition 

Some  calculations  were  performed  to  assess  the  effect  of  different  gas  compositions  on  the 
equilibrium  water  content  of  the  vapor.  This  was  done  for  each  of  the  major  components  of  the 
synthesis  gas  by  adding  0.01  to  its  mole  fraction,  and  reducing  the  other  major  components 
proportionately.  The  resulting  effects  were  quite  small.  Increasing  the  H2  mole  fraction  in  the 
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synthesis  gas  by  0.01  changed  the  water  equilibrium  mole  fraction  by  -0.00005.  Similar 
increments  for  CO  and  CO2  produced  changes  of -0.00001  and  +0.00009,  respectively.  So,  we 
can  say  that  the  water  content  in  the  vapor  is  enhanced  if  there  is  more  CO2  in  the  synthesis  gas 
and  reduced  if  there  is  more  H2  in  the  synthesis  gas,  but  these  effects  are  quite  small. 


Liquid-Phase  Composition 

While  the  water  content  of  the  equilibrated  vapor  is  the  main  quantity  of  interest  in  this  study,  it 
is  worthwhile  to  make  a few  observations  about  the  content  of  the  liquid  phase.  The  calculated 
values  for  the  different  gas  solubilities  at  the  baseline  conditions  are  shown  in  Table  2-1. 

Because  CO2  is  more  soluble  in  water  than  most  gases,  its  mole  fraction  in  the  equilibrium  liquid 
is  larger  than  those  of  the  other  gases  combined,  despite  its  being  only  18  % of  the  synthesis  gas. 
This  CO2  solubility  may  be  important  if  basic  components  such  as  NH3  are  present. 

Because  the  same  liquid-phase  fugacity  model  was  used  in  all  calculations,  the  differing  liquid 
compositions  shown  in  Table  2-1  are  attributable  entirely  to  the  effect  on  the  equilibrium  of  the 
vapor  fugacity  models.  In  particular,  the  ideal-gas  model  overestimates  the  amount  of  gas  in  the 
liquid  phase,  because  it  overestimates  the  equilibrium  mole  fractions  of  the  synthesis  gases  in  the 
vapor  phase. 

Liquid-phase  compositions  may  be  much  less  accurate  if  a different  liquid  fugacity  model  is 
used.  A few  sample  calculations  were  performed  where  the  Peng-Robinson  equation  was  used 
not  only  for  the  vapor  fugacities,  but  also  for  the  liquid  fugacities.  This  option  might  be  more 
easily  available  in,  for  example,  process  simulation  packages.  While  the  vapor-phase 
compositions  from  these  calculations  were  still  quite  accurate,  the  mole  fractions  of  the  gases  in 
the  liquid  phase  were  significantly  in  error,  by  up  to  a factor  of  two  in  some  cases.  However,  if 
the  liquid-phase  mole  fractions  are  of  no  interest,  such  an  approach  could  be  reasonable. 
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Main  Results 

Calculations  were  performed  of  the  vapor-liquid  equilibrium  between  water  and  a typical 
synthesis  gas  at  conditions  of  temperature  and  pressure  appropriate  for  a “quench”  of  the  gas.  For 
the  liquid  phase,  these  calculations  used  the  best  available  data  for  Henry’s  constants  of  gases  in 
water  and  related  quantities.  For  the  vapor  phase,  they  used  the  thermodynamically  rigorous 
virial  expansion  truncated  at  the  level  of  the  second  virial  coefficient.  When  all  uncertainties  are 
considered,  the  equilibrium  mole  fraction  of  water  in  the  vapor  phase  can  be  considered  accurate 
to  within  0.01,  with  the  largest  contribution  to  the  uncertainty  coming  from  the  neglect  of  higher- 
order  virial  coefficients. 

Other  methods  of  computing  the  vapor-phase  fugacity  were  examined.  The  ideal-gas  assumption 
seriously  underestimates  the  water  content  in  the  vapor.  The  Peng-Robinson  equation  produces 
only  a slight  underestimation,  similar  in  size  to  the  uncertainty  of  the  more  rigorous  calculation. 

It  therefore  could  be  a reasonable  option  for  engineering  calculations  for  this  system,  at  least  if 
the  form  of  the  Peng-Robinson  equation  used  is  that  documented  in  Appendix  C (in  particular, 
the  water- specific  a(T)  term,  Eq.  (C25),  is  important).  Preferably,  the  use  of  the  Peng-Robinson 
equation  in  this  context  would  only  be  for  the  vapor  phase,  with  Henry’s  law  being  used  for  the 
liquid  phase.  If  the  Peng-Robinson  equation  were  used  for  both  phases,  the  solubilities  of  gas 
components  in  the  liquid  would  not  be  accurately  represented. 

The  effects  of  temperature  and  pressure  on  the  water  content  of  the  vapor  are  as  expected;  an 
increase  of  temperature  increases  the  water  content,  while  an  increase  of  pressure  reduces  it.  The 
effect  of  gas  composition  on  the  equilibrium  water  content  is  very  small. 


Note  on  “Adiabatic”  Flash  Calculations 

While  the  scope  of  this  project  was  limited  to  flash  calculations  at  fixed  temperature  and 
pressure,  a closer  approximation  to  the  real  industrial  situation  would  be  what  is  commonly 
known  as  an  “adiabatic”  flash  calculation,  in  which  the  inlet  vapor  and  liquid  streams  are  brought 
together  at  a fixed  pressure  and  the  outlet  temperature  and  vapor  and  liquid  compositions  are 
determined  under  the  constraint  of  constant  enthalpy.  While  rigorous  adiabatic  flash  calculations 
would  be  a major  additional  effort  beyond  the  scope  of  this  work,  we  can  offer  some  qualitative 
guidance  regarding  how  the  results  of  this  work  would  translate  into  such  a situation. 


Because  the  solubility  of  the  synthesis  gas  in  water  is  small,  the  enthalpy  balance  will  be 
dominated  by  the  vaporization  of  the  water.  More  vaporization  will  result  in  more  evaporative 


3-1 


DISCUSSION 


cooling;  this  will  be  attenuated  somewhat  by  the  lower  vapor  pressure  of  the  water  as  the 
temperature  decreases.  For  a given  ratio  of  quench  water  to  synthesis  gas,  methods  that 
underestimate  the  equilibrium  water  content  of  the  vapor  at  a fixed  temperature  and  pressure 
(such  as  the  ideal-gas  assumption)  will  overestimate  the  outlet  temperature  from  the  quench  step 
in  an  adiabatic  situation. 


Possible  Enhancements 

Finally,  we  consider  some  possibilities  for  enhancement  in  the  current  calculations.  These  are 
considered  in  two  categories.  The  first  category  involves  improvements  in  the  accuracy  of  the 
calculations  described  here.  The  second  category  involves  enhancements  in  the  scope  of  the 
calculations  to  encompass  additional  aspects  not  considered  here. 

In  the  first  category,  as  mentioned  in  Chapter  2,  the  primary  source  of  uncertainty  is  the  omission 
of  higher-order  terms,  particularly  the  third  virial  coefficient  of  pure  water.  Including  these 
higher-order  contributions  could  potentially  reduce  the  uncertainty  of  the  calculated  vapor  water 
content  by  approximately  a factor  of  two.  The  most  important  task  for  such  an  enhancement 
would  be  a critical  evaluation  of  the  third  virial  coefficient  of  water,  in  order  to  estimate  its  best 
value  and  it  uncertainty.  Uncertainties  could  also  be  reduced  to  a lesser  extent  by  improvements 
in  the  theoretical  results  for  the  water-gas  second  virial  coefficients;  in  particular,  the  H2O/CO 
and  H2O/CO2  interactions  are  not  known  as  accurately  as  one  might  hope. 

In  the  second  category,  implementing  “adiabatic”  calculations  as  described  in  the  previous 
section  is  feasible,  but  would  entail  significant  additional  effort  to  superimpose  an  enthalpy 
balance  on  the  VLE  calculations.  The  inclusion  of  additional  trace  components  in  the  synthesis 
gas,  beyond  those  listed  in  Table  1-1,  could  also  be  considered.  Sufficient  data  exist  to  add  O2 
and  CH4;  parameters  for  these  gases  are  listed  in  Appendices  B and  C.  A significant  trace 
component  is  H2S.  Sufficient  data  exist  to  incorporate  H2S  in  the  liquid-phase  fugacity  model, 
but  the  H2O/H2S  interaction  for  the  vapor-phase  virial  model  is  not  known.  With  regard  to 
potential  effects  of  these  trace  components,  it  can  be  noted  that,  as  pointed  out  in  Chapter  2,  the 
water  content  of  the  equilibrium  vapor  is  not  very  sensitive  to  small  changes  in  the  synthesis  gas 
composition. 

Finally,  there  is  a significant  complication  introduced  by  the  presence  of  NH3  in  the  synthesis 
gas.  NH3  dissolved  in  the  liquid  phase  would  undergo  acid-base  reactions  with  dissolved  CO2 
and  H2S;  these  reactions  would  significantly  enhance  the  solubilities  of  these  gases.  As  a result, 
the  present  approach  would  significantly  underestimate  the  amount  of  NFI3  in  the  equilibrium 
liquid,  and  the  solubility  of  the  acid  gases  would  also  be  somewhat  underestimated.  However,  as 
long  as  the  NH3  content  in  the  synthesis  gas  was  small,  the  effect  on  the  equilibrium  water 
content  of  the  vapor  (the  main  quantity  of  interest  in  this  work)  would  be  small.  In  principle  it 
would  be  possible  to  incorporate  liquid-phase  acid-base  reactions  into  the  model,  but  it  would  be 
a major  effort. 
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Definition  of  Problem 

The  calculation  defined  for  the  project  is  the  classic  “flash”  problem,  discussed  in  many 
textbooks  on  chemical  engineering  thermodynamics  [6,7]. 

The  known  inputs  are  temperature  T,  pressure  p,  and  overall  composition  (in  this  case,  specified 
flows  of  water  in  one  input  stream  and  synthesis  gas  of  specified  composition  in  another  input 
stream).  The  output  variables  are  the  mole  fractions  of  each  component  at  equilibrium  in  both  the 
vapor  and  liquid  phases,  and  the  fraction  of  the  total  feed  in  each  phase.  The  following  equations 
must  be  solved  under  the  constraint  that  both  phases  are  at  T and  p: 

(for  each  component  0 (Al) 

z-  - ax-  + (1  - a)y.  (for  each  component  /),  (A2) 

where/  is  the  fugacity  of  component  i (which  is  a function  of  T,  p,  and  composition),  Zi,  yu  and  x, 
are  the  mole  fraction  in  the  total  feed,  the  equilibrium  vapor,  and  the  equilibrium  liquid,  and  D is 
the  fraction  (by  mole)  of  the  total  feed  in  the  liquid  phase  at  equilibrium.  Equation  (Al)  is  the 
condition  of  chemical  equilibrium  between  the  phases,  while  Eq.  (A2)  represents  conservation  of 
mass  for  each  component. 

For  phase-equilibrium  calculations,  it  can  be  more  convenient  to  work  with  the  fugacity 
coefficient,  defined  by  = fjix^p) . Equation  (Al)  then  becomes 


(for  each  component  /), 


where  the  so-called  “K-factor”  Ki  is  a key  parameter  in  the  flash  algorithm. 


(A3) 


Description  of  Algorithm 


Many  algorithms  have  been  proposed  for  solving  the  “flash”  problem;  for  systems  such  as  this 
where  the  phases  are  very  dissimilar,  a simple  algorithm  can  be  used,  consisting  of  the  following 
steps: 
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1.  Assume  an  equilibrium  composition  for  each  phase.  In  this  case,  we  initially  assume  that  the 
liquid  phase  is  pure  water,  while  the  vapor  phase  is  the  composition  of  the  synthesis  gas. 

2.  With  this  assumed  composition,  compute  the  fugacity  coefficients  and  for  each 
component  in  each  phase,  yielding  Ki  for  each  component  (see  Eq.  A3). 

3.  Perform  an  “ideal  flash,”  solving  Eqs.  (A2)  and  (A3)  with  the  values  of  Kt  fixed  at  their 
values  determined  in  Step  2.  This  is  a simple  matter  of  finding  a value  for  D that  satisfies  the 
material  balance  constraint  (A2)  and  the  constraints  that  = 1 and  = 1 . We  use  an 

i i 

ideal  flash  algorithm  originally  programmed  by  Topliss  [8]. 

4.  If  the  compositions  computed  in  Step  3 have  not  changed  (within  some  tolerance),  the 
calculation  has  converged.  Otherwise,  store  the  new  vapor  and  liquid  compositions 
determined  in  Step  3,  and  go  back  to  Step  2. 

In  the  systems  considered  here,  this  algorithm  typically  converges  within  about  five  iterations. 
We  note  that,  due  to  the  model  adopted  here  for  liquid  fugacities  (see  Appendix  B),  is 
always  independent  of  composition,  so  it  only  needs  to  be  computed  the  first  time  in  Step  2.  For 
calculations  where  an  ideal-gas  model  is  assumed  for  the  vapor,  is  also  independent  of 

composition  (because  it  is  always  1 for  all  species),  so  the  initial  ideal  flash  produces  the  final 
answer. 
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Fugacity  of  Gaseous  Solutes 

The  fundamental  thermodynamic  quantity  describing  the  solubility  of  a gas  in  a liquid  is  the 
Henry’s  constant  defined  by 


(Bl) 


where/;  and  x,  are,  respectively,  the  liquid-phase  fugacity  and  mole  fraction  of  the  solute  i.  While 
this  definition  can  be  applied  at  any  state  of  the  solvent,  in  this  work  we  restrict  our  attention  to 
vapor-liquid  equilibrium  conditions  in  single  solvents,  so  that  is  a function  only  of 
temperature  along  the  vapor-liquid  saturation  curve  of  the  solvent. 

At  small  solute  gas  pressures,  Eq.  (Bl),  with  the  limit  removed,  would  be  sufficient  to  relate  the 
solute  fugacities  to  the  mole  fraction.  However,  there  is  a (small  but  not  negligible) 
thermodynamic  correction  for  higher  pressures,  giving  the  result:  [6] 


(B2) 


where  v”  is  the  infinite-dilution  partial  molar  volume  of  solute  i in  the  solvent,  p is  the  pressure, 

is  the  saturation  vapor  pressure  of  the  solvent,  R is  the  molar  gas  constant,  and  T is  the 
absolute  temperature.  If  the  pressure  is  equal  to  the  solvent  vapor  pressure  the  correction  term 
becomes  zero  and  Eq.  (Bl)  is  recovered. 

For  phase-equilibrium  calculations  as  implemented  in  this  work,  what  is  required  is  not  the 
fugacity/;  but  the  fugacity  coefficient  = fili^iP)  • Manipulation  of  Eq.  (B2)  produces  the 

following  expression  for  the  solute  fugacity  coefficient 


(B3) 
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The  most  important  factor  in  Eq.  (B3)  is  the  Henry’s  constant  ^h-  We  use  values  from  the  critical 
data  evaluation  of  Fernandez-Prini  et  ah,  who  correlated  their  results  for  each  solute  as  a function 
of  temperature  with  a theory-based  expression:  [2] 

HK I p')=  A —f—  + C{T^ )■"“  . (B4) 

''r 


where  t=  I - Tr,  Tr  = T / Tc,  and  Tc  is  the  critical  temperature  of  the  solvent  (647.096  K for 
H2O).  Table  B-1  contains  the  parameters  A,  B,  and  C for  each  solute  in  water,  along  with  the 
minimum  and  maximum  temperatures  of  correlation. 


Table  B-1 

Parameters  for  Eq.  (B4),  along  with  Maximum  and  Minimum  Temperatures  of  Data  Used  for 
the  Correlation  [2] 


Solute 

A 

B 

C 

Tmin  (K) 

Tmax  (K) 

Ar 

-8.40954 

4.29587 

10.52779 

273.19 

568.36 

H2 

-4.73284 

6.08954 

6.06066 

273.15 

636.09 

N2 

-9.67578 

4.72162 

1 1 .70585 

278.12 

636.46 

O2 

-9.44833 

4.43822 

1 1 .42005 

274.15 

616.52 

CO 

-10.52862 

5.13259 

12.01421 

278.15 

588.67 

CO2 

-8.55445 

4.01195 

9.52345 

274.19 

642.66 

CH4 

-10.44708 

4.66491 

12.12986 

275.46 

633.11 

Both  Eqs.  (B3)  and  (B4)  require  the  vapor  pressure  of  water  at  the  temperature  of  interest. 
This  is  calculated  from  the  international  standard  correlation  of  Wagner  and  Pruss,  which  has  the 
following  functional  form: 

ln(/>"  / (^c  ! T){a^T  a2T^'^  +a^T^'^  (B5) 

where pc  is  water’s  critical  pressure  (22.064  MPa),  and  rand  Tc  are  defined  below  Eq.  (B3).  The 
values  of  the  coefficients  are  ai  = -7.85951783,  <32  = 1.84408259,  <23  = -11.7866497, 

(34  = 22.6807411,  as  = -15.9618719,  = 1.80122502.  [91 

The  final  quantity  required  for  computation  of  the  solute  fugacity  coefficients  by  Eq.  (B3)  is  the 
solute  partial  molar  volume  i7”  at  infinite  dilution  in  water.  This  quantity  is  calculated  as 

described  by  Fernandez-Prini  and  Japas  based  on  a molecular  perturbation  theory.  [10,11]  Their 
final  expression  is 


vr-RTzM  + 


\-T] 


1 + 3 


r + r 


+ 9- 


Tjr^ 


\-r]  {\-T])  {\-p) 


(B6) 
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where  D is  the  “packing  fraction”  of  the  pure  solvent  (water),  defined  by  ;;  = TiN^p^dl  / 6 , 

where  is  the  molar  density  of  water,  is  the  characteristic  diameter  of  a water  molecule 
(taken  as  0.27  nm  in  all  cases  here),  and  Na  is  Avogadro’s  number,  r is  the  ratio  where  <7,  is 
a characteristic  diameter  for  the  solute  molecule.  Dr  is  the  isothermal  compressibility  of  pure 
water. 


The  pure-water  properties  p^  and  Dr  are  obtained  from  the  NIST  Standard  Reference  Database 
for  water  properties.  [12]  The  only  parameter  required  for  the  calculation  of  v”  is  then  which 
is  correlated  as  a function  of  temperature  by  one  of  two  functional  forms:  [11] 


d,  =d°  +Z),.(7’-298K), 

(B7a) 

d.  =d°  +Z);in(7’/298K). 

(B7b) 

Table  B-2  contains  the  parameters  d°  and  bi  or  b[  for  all  the  solutes  considered  in  this  work. 


Table  B-2 

Parameters  for  Eq.  (B7)  (last  column  is  bi  for  all  solutes  except  CO2,  where  it  is  b[ ) [11]. 


Solute 

d°  (nm) 

bi  (nm/K)  or  b'  (nm) 

Ar 

0.241 

-2.01  X 10^ 

H2 

0.266 

-1.25  X 10'^ 

N2 

0.339 

-1.70  X 10"^ 

O2 

0.342 

-2.00  X 10^ 

CO 

0.353 

-2.53  X 10'^ 

CO2 

0.427 

-0.2246 

X 

0 

0.369 

-3.05  X 10'^ 

We  note  that,  at  the  pressures  and  temperatures  of  interest  in  this  work,  the  correction  to  the 
solute  fugacity  coefficients  provided  by  the  v ” term  is  typically  on  the  order  of  one  or  two 
percent.  Its  accuracy  is  therefore  not  crucial  in  the  context  of  the  present  project. 

Fugacity  of  Liquid  Water 

Under  the  conditions  considered  here  (dilute  solution  of  solute  gases),  the  calculation  of  the 
fugacity  of  liquid  water  is  straightforward,  given  merely  by  its  fugacity  in  the  pure  state  at 

the  temperature  and  pressure  of  interest  (obtained  from  [12])  multiplied  by  its  mole  fraction 
(which  is  slightly  less  than  1 due  to  the  solubility  of  the  solute  gases).  The  equation  for  the 
fugacity  coefficient  becomes 
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P 


Impact  on  Calculations 

It  can  be  noted  that  the  main  quantity  of  interest  in  this  project  (the  amount  of  water  in  the 
equilibrium  vapor  phase)  is  not  very  sensitive  to  the  liquid-phase  fugacity  calculations  (other 
than  the  calculation  of  , which  is  well-known).  The  solute  fugacity  coefficients  as 

calculated  in  Eq.  (B3)  primarily  influence  the  solubility  of  the  gases  Xi.  While  these  solubilities 
may  be  of  interest  for  some  applications,  they  are  small  enough  that  their  impact  on  the 
equilibrium  fugacity  of  water  is  generally  less  than  1 %. 

We  also  note  that,  due  to  the  structure  of  the  model,  both  the  solute  fugacity  coefficients  (Eq. 
B3)  and  the  water  fugacity  coefficient  ^ (Eq.  B8)  depend  only  on  temperature  and  pressure,  but 
not  on  the  liquid-phase  composition.  This  simplifies  the  vapor-liquid  equilibrium  calculations,  as 
discussed  in  Appendix  A. 
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Calculation  from  Virial  Coefficients 

For  gas-phase  systems,  a rigorous  description  of  the  thermodynamics  is  provided  by  the  virial 
expansion,  which  is  written  as  a series  of  corrections  to  the  ideal-gas  law: 


— ^ -\  + Bp  + Cp^  + . . . , 

pRT 


(Cl) 


where  p is  the  pressure,  p the  molar  density,  R the  molar  gas  constant,  and  T the  absolute 
temperature.  The  quotient  z (which  is  1 for  an  ideal  gas)  is  usually  called  the  compressibility 
factor,  B is  called  the  second  virial  coefficient,  C the  third  virial  coefficient,  and  so  forth.  An 
alternative  form  of  the  expansion  can  be  formulated  as  a power  series  in  p rather  than  p,  but  it 
has  been  shown  that  the  pressure  form  is  generally  less  accurate  than  the  density  form  when 
truncated  at  the  same  level.  [13] 

For  mixtures,  the  virial  coefficients  are  rigorously  given  as  mole-fraction-weighted  sums  over 
characteristic  binary  quantities  (for  5),  ternary  quantities  (for  C),  etc.  For  B,  the  expression  is 


(C2) 


where  y is  the  mole  fraction  of  the  species  in  the  vapor  mixture  and  the  sums  extend  over  all 
components.  Bij  is  completely  determined  by  the  interactions  between  one  molecule  of  species  i 
and  one  of  species  j.  For  a given  pair,  it  is  a function  of  temperature  only.  If  iff,  BifT)  is  just  the 
pure-component  value  of  5(7);  these  values  are  well  known  for  the  substances  of  interest  here. 
Obtaining  reliable  values  of  BifT)  for  unlike  pairs  is  more  difficult,  and  has  been  the  main  focus 
of  this  project. 

At  moderate  pressures,  it  is  adequate  to  truncate  Eq.  (Cl)  after  the  term  involving  the  second 
virial  coefficient.  Standard  thermodynamic  manipulations  then  produce  the  following  expression 
for  the  fugacity  coefficient  for  each  species:  [6] 


(C3) 


where  the  sum  is  over  all  species  including  i.  The  molar  density  p (and  therefore  the 
compressibility  factor  z)  is  solved  for  iteratively  from  Eq.  (Cl)  once  all  the  virial  coefficients  are 
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known.  Use  of  Eq.  (C3)  requires  knowledge  of  all  the  second  virial  coefficients  for  all 
pairs  of  species  in  the  system. 

B(T)  Data  for  Pure  Components 

B{T)  for  pure  water  is  taken  from  the  correlation  of  Harvey  and  Lemmon,  which  has  the  form 
B{T)IB^  = a, +a^{ry^-^^  +a^{r)-^\  (C4) 

where  5°  = 1000  cmVmol,  T = 7/100  K,  a\  = 0.34404,  a2  = -0.75826,  a^,  = -24.219,  and 
^4  = -3978.2.  [14] 

B{T)  for  pure  oxygen  is  taken  the  work  of  from  Wagner  et  al.: 

5(7)75"  =/7,  +b^{ry^-^^  +/?3(7*)-''  +b,{ry^\  (C5) 

where  5°  = 1000  cmVmol,  f = Til  K,  by  = 0.143389,  b2  = -0.629863,  b^  = -0.577814  x 10^ 

Z?4  = 0.695858  x 10^  and  ^5  = -0.246023  x lo".  [15] 

5(7)  for  pure  hydrogen  is  taken  from  a correlation  by  Hodges  et  al.: 

5(7)/ 5"  =ai(7*)'"”  + a^yr*)-'-^  + a^{ry^-^  + a^{ry^\  (C6) 

where  5^  = 1 cmVmol,  y = 7/100  K,  a\  = 42.0803,  a2  = -143.982,  a^,  = 146.918,  and 
^4  = -47.5601.  [16] 

For  other  pure  species,  5(7)  is  taken  from  the  reference-quality  equations  of  state  implemented 
in  the  NIST  Refprop  database.  [3]  The  expressions  for  5(7)  derived  from  these  equations  take 
the  form: 

S(7’)/A=S«<^'’''  <C7) 

i 

where  is  the  critical  density  of  the  fluid,  r = TJT,  where  7c  is  the  critical  temperature  of  the 
fluid,  and  the  a,  and  by  are  a series  of  parameters  specific  to  the  fluid.  Table  C-1  gives  parameters 
fl,  and  bi  for  use  in  Eq.  {Cl),  and  Table  C-2  lists  the  values  of  the  critical  constants. 
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Table  C-1 

Parameters  for  Eq.  (C4)  for  S(7)  of  Ar,  N2,  CO,  CO2,  and  CH4 


Param. 

Ar 

N2 

CO 

0 

0 

X 

0 

a^ 

0.08872230499 

0.924803575275 

0.90554 

0.388568232032 

0.04367901028 

32 

0.705148051673 

-0.492448489428 

-2.4515 

2.93854759427 

0.6709236199 

-1.68201156541 

0.564857472498 

0.53149 

-5.58671885349 

-1 .765577859 

34 

-0.149090144315 

-1.61720005987 

-0.12778 

-0.767531995925 

0.1830487909 

as 

-0.120248046009 

-0.481395031883 

(not  used) 

2.16589615432 

0.1511883679 

as 

-0.321813917507 

-0.0435762336045 

(not  used) 

-0.370556852701 

-0.4289363877 

aj 

0.332300176958 

(not  used) 

(not  used) 

-0.0167758797004 

-0.01932040831 

hi 

0 

0.25 

0.25 

0 

-0.5 

bz 

0.25 

0.875 

1.125 

0.75 

0.5 

bs 

1 

0.5 

1.5 

1 

1 

b4 

2.75 

0.75 

3.625 

2 

0 

bs 

4 

2 

(not  used) 

1.5 

1 

be 

3 

4 

(not  used) 

3 

2 

by 

3.5 

(not  used) 

(not  used) 

6 

5 

Table  C-2 

Critical  Parameters  and  Acentric  Factor  for  Substances  Considered  in  this  Work 


Substance 

7‘c(K) 

Pc  (kPa) 

Pc  (mol/dm^) 

0) 

H2O 

647.096 

22064. 

17.873728 

0.3443 

Ar 

150.687 

4863. 

13.4074 

-0.00219 

H2 

33.145 

1296.4 

15.508 

-0.219 

N2 

126.192 

3395.8 

11.1839 

0.0372 

O2 

154.581 

5043. 

13.63 

0.0222 

CO 

132.86 

3494. 

10.85 

0.0497 

CM 

0 

0 

304.1282 

7377.3 

10.6249 

0.22394 

CH4 

190.564 

4599.2 

10.139 

0.01142 

B(T)  Data  for  Aqueous  Binaries 

The  for  pairs  where  water  is  one  of  the  components  are  the  most  important  factors 
determining  the  vapor-phase  fugacity  coefficient  in  Eq.  (C3),  and  therefore  are  the  most 
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important  in  determining  the  equilibrium  mole  fraction  of  water  in  the  vapor.  For  the  most  part, 
these  are  obtained  from  a series  of  recently  published  works  making  use  of  quantum  chemistry. 
[16,17,18,19,20,21]  The  exception  is  the  H2O/CO2  system,  where  the  experimental  data  of  Patel 
et  al.  cover  the  temperature  range  of  interest  and  were  fitted  for  this  work  to  a simple  functional 
form.  [4] 

The  equations  for  Bij(T)  for  most  of  these  pairs  take  a form  similar  to  that  of  Eq.  (C4): 

^ a, (T*)^  +a,{rf^  +a,{rf\  (C8) 

where  B^  = 1000  cmVmol  and  T*  = 77100  K.  Values  of  the  a and  b coefficients  for  each  aqueous 
binary  for  use  in  Eq.  (C8)  are  given  in  Table  C-3.  Information  about  the  uncertainty  of  these 
correlations  may  be  found  in  the  original  references. 


Table  C-3 

Parameters  for  Eq.  (C8)  for  BifJ)  for  H2O  with  Various  Gases 


Parameter 

HaO/Ar  [17] 

H2O/H2  [16] 

H2O/N2  [19] 

H2O/O2  [20] 

H2O/CO  [21] 

H2O/CO2  [4] 

ai 

96.1941 

33.047 

67.595 

124.605 

485.062 

-3308.1 

32 

-211.074 

-250.41 

-249.83 

-214.421 

-2054.24 

(not  used) 

^3 

-96.4425 

285.42 

-204.38 

-102.818 

7422.08 

(not  used) 

34 

-12.6006 

-186.78 

(not  used) 

-22.360 

-5919.89 

(not  used) 

bi 

-0.31 

-0.21 

-0.24 

-0.33 

-0.18 

-2.66 

bz 

-0.82 

-1.50 

-1.06 

-0.73 

-1.70 

(not  used) 

63 

-2.24 

-2.26 

-3.22 

-2.03 

-2.10 

(not  used) 

64 

-4.60 

-3.21 

(not  used) 

-4.07 

-2.25 

(not  used) 

The  only  aqueous  binary  for  which  5y(7)  is  not  described  by  Eq.  (C8)  is  H2O/CH4,  for  which  the 
correlation  is 

fi„(7’)/(l  cm’  /mol)  = f^aXTjT)  , (C9) 

/=1 

where  Tq  = 599.216  K,  ai  = 109.447,  a2  = -221 .69 A,  = 182.497,  = -80.944,  as  = 18.368, 

and  «6  = -1.673.  [18] 

B(T)  Data  for  Nonaqueous  Binaries 

The  Bij(T)  for  pairs  where  neither  component  is  water  play  a minor  but  not  negligible  role  in  the 
phase  equilibrium.  For  these  binaries,  we  make  use  of  a corresponding-states  correlation 
developed  by  Tsonopoulos.  [22]  The  correlation  takes  the  form 
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B{T)p^ 

RT^ 


_ ^(0) 


T 

V^c 


+ 0)F 


(1) 


r j,\ 


kT.j 


(CIO) 


with 


’(0) 


f jS 


yT.j 


= 0.1455 


0.330  0.1385  0.0121  0.000607 


’(1) 


yT.j 


= 0.0637  + 


0.331  0.423  0.008 


(Cll) 

(C12) 


Pc  and  Tc  are  the  critical  pressure  and  temperature,  respectively,  co  is  the  acentric  factor,  a semi- 
theoretical  parameter  related  to  the  degree  to  which  a substance  deviates  from  a “simple”  fluid; 
values  for  co  are  well-known  for  most  common  fluids.  Tr  is  the  reduced  temperature,  7r  = TUc. 
Critical  properties  and  acentric  factors  for  all  fluids  used  in  this  work  are  taken  from  the  NIST 
Refprop  database  [3]  and  are  listed  in  Table  C-2. 

For  “quantum”  gases  (those  such  as  H2  with  such  low  molecular  masses  that  their 
thermodynamic  properties  are  affected  by  quantum  mechanics),  temperature-dependent 
“effective”  critical  constants  are  used,  defined  as  follows: 


Z = 


1 + 


21.8  ’ 
MT 


Po  = 


Pc 


1 + 


44.2 

MT 


(C13) 


(C14) 


and  pi  are  hypothetical  values  that  the  critical  temperature  and  pressure  would  have  in  the 
absence  of  quantum  effects,  and  M is  the  relative  molar  mass.  For  H2,  the  only  quantum  fluid 
considered  in  this  project,  = 43.6  K,  pi  = 2.05  MPa,  and  M = 2.01594. 


In  order  to  use  Eq.  (CIO)  to  derive  Bij{T)  for  unlike  pairs,  combining  rules  are  required  to  obtain 
characteristic  values  of  pc,  co,  and  M (if  one  of  the  pair  is  H2)  for  use  in  the  Tsonopoulos 
correlation.  As  suggested  in  the  literature,  the  following  combining  rules  are  used:  [6,22] 

(C15) 


Pc,tJ  = 


ITc,i+ Pcj^cj  IT,j) 


(C16) 
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cOy  =0.5[co^+o)j). 

(C17) 

My'  -0.5(m:’  +m;‘). 

(C18) 

The  binary  parameter  kij  in  Eq.  (C15)  is  set  to  zero  in  the  absence  of  data,  but  can  be  adjusted  to 
fit  experimental  data  if  they  are  available.  Table  C-4  lists  values  of  kij  (which  is  the  same  as  %) 
and  the  corresponding  sources  of  mixture  data  for  a number  of  binaries  in  this  work.  For  those 
not  listed  in  the  table,  a value  of  0 is  used. 


Table  C-4 

Binary  Parameters  kij  ior  Eq.  (C15)  for  Pairs  of  Nonaqueous  Gases 


Component  / 

Component  j 

f<ij 

Data  Source 

Ar 

H2 

-0.12 

[23] 

Ar 

N2 

-0.01 

[23] 

Ar 

O2 

0.00 

[23] 

Ar 

CO 

-0.01 

[24] 

Ar 

CM 

0 

0 

0.05 

[23] 

Ar 

CH4 

0.02 

[23] 

H2 

N2 

-0.04 

[24] 

H2 

CO 

-0.04 

[24] 

H2 

C02 

-0.04 

[25] 

H2 

CH4 

-0.05 

[24] 

N2 

02 

-0.01 

[23] 

N2 

CO 

-0.02 

[26] 

N2 

C02 

-0.01 

[23,27] 

N2 

CH4 

0.03 

[23,28] 

O2 

C02 

0.05 

[23] 

O2 

CH4 

0.04 

[23] 

CO 

0 

0 

-0.04 

[29] 

CO 

CH4 

0.01 

[24] 

CM 

0 

0 

CH4 

0.05 

[23,27] 
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Supplementary  Calculation  with  Peng-Robinson  Equation  of  State 

For  comparison  purposes,  some  calculations  for  this  project  were  performed  with  the 
semiempirical  Peng-Robinson  equation  of  state.  [30,311  The  Peng-Robinson  equation  has  the 
following  form; 


RT a{T) 

v-b  v(v -\-b)  + b{v - b)  ' 


(C19) 


where  p is  the  pressure,  T is  the  absolute  temperature,  v is  the  molar  volume,  and  R is  the  molar 
gas  constant,  a and  b are  substance-specific  parameters;  Z?  is  a constant  depending  only  on  the 
critical  properties,  while  a is  a function  of  temperature  in  a manner  that  depends  on  the  acentric 
factor  co: 


RT 

6 = 0.778 — ^ 


a{T)  = a{TJa{T), 


«(?;)  = 0.45724 


{RTy 


«(r)  = [l  + /?(l-V7V7^f , 


p = 0.37464  + 1 .5422661  - 0.269926i' 


(C20) 

(C21) 

(C22) 

(C23) 

(C24) 


We  make  use  of  a special  parameterization  introduced  by  Peng  and  Robinson  for  a{T)  for  water: 

[31] 


{T)  = [l .0085677  + 0.821 54(l  - / 7;  f , for  ^T/T;  <0.85. 


(C25) 


When  yjT  IT^  > 0.85 , Eqs.  (C23-C24)  are  used  as  with  other  fluids. 

Use  of  the  Peng-Robinson  equation  for  mixtures  requires  mixing  rules  for  the  parameters  a and  b 
based  on  their  pure  component  values: 

« = , (C26) 

> i 

b = , (C27) 
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where  yt  is  the  mole  fraction  of  component  i.  For  unlike  pairs,  a combining  rule  is  needed  for  aij 
based  on  pure-component  values  a,  and  af. 

a,i=[a,ay\  (C28) 


It  is  possible  to  introduce  an  adjustable  binary  parameter  into  Eq.  (C28)  (similar  to  that  in  Eq. 
(C15))  in  order  to  attempt  to  fit  experimental  data  more  closely,  but  that  enhancement  was  not 
explored  in  this  project. 

From  the  Peng-Robinson  equation,  the  equation  for  the  fugacity  coefficient  of  a component  i 
is:  [6] 


pv 

RT 


-1 


-In 


a 


P(v-b) 

RT  l^lbRT 


a 


h 

b 


In 


(C29) 


The  accuracy  of  our  implementation  of  the  Peng-Robinson  equation  was  tested  by  verifying  that 
calculated  results  for  densities  and  fugacity  coefficients  matched  those  from  an  independently 
programmed  implementation  available  in  a NIST  database.  [3] 
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